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Abstract 

Floquet theory, first published in 1883 for periodic linear differen¬ 
tial equations, is extended in this paper to multitime diagonal recur¬ 
rences. We find explicitly a monodromy matrix, and we comment its 
eigenvalues (called Floquet multipliers). The Floquet point of view 
brings about an important simplification: the initial linear diagonal 
recurrence system is reduced to a linear recurrence system, with con¬ 
stant coefficients along “diagonal lines”. The results are applied to 
the discrete multitime Samuelson-Hicks models, with constant, respec¬ 
tively multi-periodic, coefficients, in order to find bivariate sequences 
with economic meaning. For constant coefficients case, it was found 
also the generating function; for multi-periodic coefficients case, we 
determined the Floquet multipliers. 
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1 Discrete multitime recurrences 

The multivariate recurrences are based on multiple sequences and come from 
areas like analysis of algorithms, computational biology, information theory, 
queueing theory, hlters theory, statistical physics etc. 

We consider the lattice of points with integer positive coordinates in M”. 
A multi-variate recurrence is a set of rules which transfer a point into another, 
together with initial conditions, capable to cover the hole lattice. 

A way to study a multiple sequence is to apply a transform (generating 
function, Mellin transform etc), yielding a multi-variable complex function 
F. We then study the singularities of A, its Taylor development around 
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origin, and extract information on the original multiple sequence via multi- 
variable complex analysis (especially the explicit formula of the multivariate 
sequence). The advantages of generating functions (GFs) are: provide a 
compact representation of the multivariate sequence, are well adapted for 
computer algebra manipulation (symbolic computation avoids numerical er¬ 
rors and can yield insight); can prove identities and simplify sums; can solve 
recurrences (difference equations); can yield recurrences, and hence fast com¬ 
putation algorithms; allow extraction of asymptotic approximations. The 
multivariate GF is the best all-round tool for studying multi-variate recur¬ 
rences. 

Bousquet-Melou and Petkovsek [Tj showed that while in the univariate 
case solutions of linear recurrences with constant coefficients have rational 
generating functions (see also 0), the multivariate case is much richer, even 
though initial conditions have rational generating functions the correspond¬ 
ing solutions can have generating functions which are not rational, algebraic, 
or even D-£nite GF in two or more variables. 

Pemantle and Wilson [5] describe a symbolic method for computing gen¬ 
erating functions for multivariate recurrences. 

Hauser and Koutschan [5] give theorems of division for formal power se¬ 
ries, which clarihes the structural background and provides short, conceptual 
proofs for generating functions. In addition, extending the division to the 
context of differential operators, the case of recurrences with polynomial co¬ 
efficients can be treated in an analogous way. 

Floquet theory, hrst formulated for periodic linear ODEs [3], was ex¬ 
tended to difference equations (13, PI, i, 0). We have extended this 
theory to the multitemporal first order PDEs PI and now to multivariate 
multi-periodic recurrences, borrowing mathematical ingredients from some 
papers ra-pi involving multitime (multivariate time) dynamical systems. 
In Floquet theory it is necessary to hnd explicitly the associated monodromy 
matrix and its eigenvalues (called Floquet multipliers). 

2 Linear discrete multitime diagonal 
recurrence with periodic coefficients 

An element t = (f^,..., P") G N™' is called discrete multitime. A function of 
the type x : N™ M" is called multivariate sequence. Also, for convenience, 
we denote p(t) = min{f^, ..., t™} and 1 = (1,1,.. ., 1) G M™'. 

Let m>2 and A: —)■ AT„(M). Our aim is to continue the study of a 
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linear discrete multitime diagonal recurrence system (see 0) 


x{t + 1) = A{t)x{t), \/t G 


( 1 ) 


In the paper [1] one proves the next result. 

Theorem 1. Let m > 2, A: —>■ b: hJ™ M"'. We consider the 

(m — 1)-sequences /i, / 2 , • ■ ■, /m: ^ M"", such that 




t“=0 


( 2 ) 


W,..., r-\ r +\..., tf^-\ t ^+\..., t’" e N, 

for any a, jd E {1,2, Then the unique m-sequence x: N’" —)■ 

which verifies 

x{t + 1) = A{t)x{t) + 6(t), Vt G N™', 

ti^=o 

V/3 G {1,2, 

is defined either by the formula 


x(t) = A{t — l)A{t — 2 ■ 1) ■ ... • A{t — t^ ■ 1)- 
■ fy{A -tf^,...Sh,..., t^-^ - + b{t - 1) 

tP 

+ ^ 2l(t - l)A{t -2-1) ■ A{t-{k-l) ■ l)b{t -k-1), 

k=2 

if fi{t) =t^ >2, 


or by the formula 

x{t) = A{A - 1,, t^-^ - 1, 0, - 1,... - !)■ 

■ fy{t^ - l,...,t^,...,t'”“^ - 1) 

+ b{t^ - 1 ,..., t^~^ - 1 , 0 , - 1 ,... - 1 ), 

if fi{t) =t^ = 1. 


□ 


The function 




—A:-l), if /i(f) > 1, 

k=l 

In, if ti{t) = 0, 
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is called the fundamental matrix {transfer matrix) associated to the recur¬ 
rence ([T]). In fact, <!>(■) is the unique matrix function which verihes the 
problem 


d)(t + 1) = A(t)<h(f), WteN^, 

$(f) ^ =4, y{t\...,td,---,n^^^-\ 

t/3=0 

V/3 e {1,2,.. .,m}. 


(3) 


This follows applying the Theorem [T] for n (vector) recurrences to whom is 
equivalent the problem (|3]) (one applies the Theorem [1] for each column of 
the matrix $(■)). 

If the functions 

verify the relations ([2]), then one observes (Theorem [1]) that the unique m- 
sequence x : N™ —?■ M"" which verihes the recurrence ([1]) and 




m—1 


V/3 e (1,2,. 


can be written in the form 

x{t) = <h(f)4(t^ - ..., 4,..., - t^), if n{t) = 


If the matrix function A(-) is constant, i.e., A{t) = A, Vf, then the fun¬ 
damental matrix becomes <I>(f) = A^^^\ 

A linear discrete multitime diagonal recurrence is called T-diagonal-periodic 
{T e N*) if 

A{t + T-1) = A{t), Wt (4) 

This is the only multi-periodicity compatible to the diagonal recurrence, in¬ 
dependently of the initial conditions. 

Proposition 1. Suppose the recurrence o is T-diagonal-periodic and we 
introduce the function 

T-l 

A: A{t) = Y[A{t+{T-l-k)-l), Vt G M™. 

k=0 

Then, for any k E N*, the fundamental matrix satisfies 

$(f + fcT-l) =<!>{t)(^A{t- ^{t)-l)y. (5) 
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Proof. If = n{t), then 


j^kT = min{t^ + kT, + kT,... + kT}. 

Induction after k. For k = 1: 

The case = 0: 

+ T ■ l) = /l(t + (T - 1) ■ l) A(t + (T - 2) ■ l) ■... • A(t + (T - T) ■ l) 

= 4l(t) = Ht)A{t). 

Case > 1: 

^(t + T-l) = A(t + T-l-l)A(t+ T-1-2-1) - .•A(t + r-l-(t^ + T)-l) 

= + T • 1 - l)A{t + T-l-2-l)-...-A(t + T-l-t^-l)- 

■A{t-t^■l + {T-l)-l)A{t-t^■l + {T-2)-l)-.. .-A^t-t^■l + l)A{t-t^-l) 
= A{t - l)A{t -2-l)-...-A{t-t^-l)- 
■A{t-t^■l + {T-l)-l)A{t-t^■l + {T-2)-l)-.. .■A{t-t^■l + l)A{t-t^-l) 

= ^{t)A{t-tf^ ■ 1 ). 

Suppose the relation ([S]) is true for any k, and we shall verify for k + 1: 
$(t+(fc+l)r-l) = $(t+fcT-l+T-l) = ^t+kT-l)A{t+kT-l-{t^+kT)-l) 

/ \ Aj / \ Ai I ^ 

= <^{t)(^A{t - ■ l)] A{t-t^ -l) =$(t)M(t-t^-l)j . 

□ 

Suppose that we are in the conditions of the Proposition [TJ The matrix 
function 


D; A1„(R), D{t) = A{t- Vt G 

is called the monodromy matrix associated to the T-diagonal-periodic recur¬ 
rence For fc = 1, the formula (l5|) can be written 

+ T-1) = <I>(t)D(t), VteN’". (6) 

We observe that the relation D(t + 1^ = D(t), Vf G N™' holds. 

Moreover, let us suppose that, for any t G N™', the matrix A{t) is invert¬ 
ible, hence A{t) is also invertible; it follows that for each t G M"*, there exists 
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B{t) G Ain{C) (which is not unique), such that = A{t). For each t, 

we hx such a matrix B{t)] obviously the matrix B{t) is invertible. 

Dehne the function 

B{t) = B{t- fi{t) -1), WteN^. (7) 

We observe that the relation 

B{t + l)=B{t), VteM’” (8) 

is true, i.e., the matrix B verihes a special recurrence. It follows immediately 
the relation 

D{t) = B{tf, yt e N™. (9) 

Moreover, one observes that if for any t, the matrix A{t) is invertible, then 
the matrices <h(t) and D{t) are invertible too. 

Theorem 2. Let m > 2 and A\ N™' —with A{t) invertible, Vf G 
M™'. Suppose there exists an integer T > 1 such that the relation (jl]) is true. 

Then there exists a function P: N™ —)■ Ain{C), with the property that 
P{t + T-1) = P{t), yt, and such that the fundamental matrix of the T-multi- 
periodic recurrence is written 

^t) = P{t)B{tY^^\ VtGN’" (10) 

{^where B{-) is the function definite by the formula ([7])) . 

Proof. Let 

P, ^Mn{C), 

Pit) = $(t)5(f )-™°^*'yt G N’". 

It is sufficient to show P(t + T ■ 1) = Pit), Vf G 

Let t = iA,P,, t^) G N™ and t^ = /i(f). Obviously, 

t^ + T = min{f^ + T, + T,..., r + T}, 

P(t + T ■ 1) = $(f + T ■ l)Bit + T ■ 1)-^^-^. 

From the relation ([8j), we deduce Bit + k-1) = Bit),yk G N, and particularly 
Bit + T ■ 1) = Bit). From the relations ([6]), ([^, it follows <I>(f + T ■ 1) = 
<h(f)i?(f)^. Hence 

Pit + T-1) = ^t)Bit)'^Bit)-^^-'^ = ^t)Bit)-^^ = Pit). 


□ 
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In the condition of the Theorem [21 from the formula fITUD and the fact 
that the matrix <h(f) is invertible, it follows that P{t) is also invertible. 

The most important result of Floquet type is 

Theorem 3. Let m > 2 and A: —)■ with A{t) invertible, Vf G 

N™'. Suppose that there exists an integer T > 1 such that the relation (|1]) is 
true. Let B : N™' —)■ A^„(C) be the function definite by the formula ([7]). 

We consider the recurrences 

x{t + 1) = A{t)x{t), Vf G N™'; (11) 

y{t+l) = Bit)y{t), yteN^. (12) 

If y{t) is a solution of the recurrence flT2|) . then x{t) := P{t)y{t) is a 
solution of the recurrence fITT]) . And conversely, if x{t) is a solution of the 
recurrence dm), then y{t) := P{t) ^x{f) is a solution of the recurrence flT2]) . 

Proof. Let y{t) be a solution of the recurrence flT2|) and x{t) := P{t)y{t); 
hence y{t) := P{t)~^x{t). We £x f G M”* and t^ = p[f). It follows 

y{t + 1) = B{t)y{t) P{t + l)“^x(f + 1) = B{t)P{t)~^x{t) 

x{t + 1) = P{t + l)B{t)P{t)~^x(t). 

We use the formula (ITU]) and deduce that the foregoing relations are equiva¬ 
lent to 


x{t + 1) = <!>{t + l)B{t + l)-^^-^B{t)B{ty^<I>{t)-^x{t). 

This relation is equivalent to (according the formula flHl) ) 

x{t -I- 1) = ^(t -I- l)B{t)~^^~^B{t)B{tY'^^{t)~^x{t) 

> x{t -I- 1) = A{t)^{t)^{t)~^x{t) 
A{t)x{t). 

□ 

In this way, the Floquet point of view brings about an important simpli- 
hcation: the linear diagonal recurrence system ca is reduced to the linear 
diagonal recurrence system flT2|) . with constant coefficients along “diagonal 
lines”. 


x{t -I-1) = <F(t -I- l)*h(f) ^x{t) 

x{t 1) = 

Like it proves the reverse. 
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3 Two-time recurrence of way required 

Let t G A first order linear bi-recurrence of the form 

x{t^ + 1 , 0 ) = Fi((t\ 0 ),x(f\ 0 )), + 1 ) = 

with initial condition a;( 0 , 0 ) = xq, is called two-time recurrence of way re¬ 
quired. Of course, similar to this recurrence, we can consider discrete multi¬ 
time recurrences of way required or discrete multitime nonholonomic recur¬ 
rences. 

Suppose we have linear two-time recurrence of way required, with con¬ 
stant coefficients, 

x{t^ - 1 - 1 , 0 )= Aix{t^, 0 ), x{A,f^ + 1 ) = A2x{t^,t‘^)- 

Theorem 4. The solution (imposed by the path) of a linear two-time recur¬ 
rence of way required, with constant coefficients, is 

x{A,f^) = A2 A-f xq. 


4 Discrete multitime constant coefficients 
Samuelson-Hicks model 


We assume that t = {A,..., f^) G N™ is a discrete multitime. Having in mind 
the discrete single-time Samuelson-Hicks model [S], we introduce a discrete 
multitime Samuelson-Hicks like model based on the following ingredients: 
(i) two parameters, the hrst 7 , called the marginal propensity to consume, 
subject to 0 < 7 < 1 , and the second a as decelerator if 0 < a < 1 , keeper 
if a = 1 or accelerator if a > 1 ; (ii) the multiple sequence Y{t) means the 
national income and is the main endogenous variable, the multiple sequence 
C{t) is the consumption] (iii) we assume that multiple sequences Y{t), C(t) 
are non-negative. 

We propose a first order discrete multitime constant coefficients Samuelson- 
Hicks model either as hrst order diagonal recurrence system 


/K(t + iA 

\C{t + 1)/ 


^7 + a 

V 7 



t G N™, 




= Yoo.{t\ ..., ..., n, C{t)\tc^o = Coa{t\ ..., ..., n 

Yo^{t\ ...,P,..., = yoy{t\ ..., 9 ,..., r)|i.=o 
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(13) 

(14) 


or as the second order homogeneous diagonal recurrence equation 

Y{t + 2) - (7 + a)Y{t + 1) + aY{t) = 0, t G N™, (16) 

y(t)|i.=o = h))a(t\...,t“,...,n, r(t)|t.=i = Fi„(t\...,P,...,n (17) 

Y^^{t\ tY n\t^=o = Yop{t\ tP, (18) 
yi„(t\P,r)|*,=i = Y^p{t\ r)|i.=i (i9) 

yoa(^^...,^“,...,nl^/^=l = n/3(ts...,t^,...,nit“=o. (20) 

4.1 Samuelson-Hicks discrete multitime recurrence 
equation with constant coefficients 

4.1.1 Characteristic equation 

Let us consider the second order homogeneous discrete multitime diagonal 
recurrence equation fllOp . with constant coefficients 7 and a. Looking for a 
solution of the form 


y(t) = < e,l >= 1, 

it appears the characteristic equation 

— (7 y o)A 'L o = 0. 

The discriminant of this equation is 5 = (7 + a)^ — la. If 5 > 0, then we 
have real roots Ai > 0, A 2 > 0; if 5 = 0, then Ai = A 2 > 0; if h < 0, then the 
roots Ai, A 2 are complex conjugated with modulus \foL. The stability requires 
|A,|<1, z = l,2. 

4.1.2 Generating function 

First computation variant We start with the recurrence 

lrn+2n+2 (d T n_|_i T cYY^an 0) ^ ^ 0, ?7. ^ 0. 

Our aim is to determine a closed form for the corresponding generating func¬ 
tion 

F{x,y)= Y. YranX^y- 

m=0,n=0 
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(analytic in a neighborhood of the origin; it can often be extended by ana¬ 
lytic continuation) in terms of the given initial conditions and the recurrence 
relation, and then to derive a direct access formula for the coefficients Ymn 
of the development of F{x,y) in Taylor series around origin. 

Multiplying with and summing, it follows 

m,n>0 m,n>0 

+a = 0. 

m,n>0 

We replace each double inhnite sum in this equation by an algebraic expres¬ 
sion involving the generating function, respectively 

y^+2n+2X”^+22/"+2 = F{x, y)-Yl ^-0^" - E ^0-^" + ^00 

m,n>0 m>0 n>0 

-y ^ ^ X] + ^y^n, 

m>l n>l 

= xy f F{x, y) + Y,o~Y. ^^nox^ “ E 

m,n>0 \ m>0 n>0 

Ymnx"'^‘^y'^^‘^ = x^y‘^F{x,y). 

m,n>0 

We denote 

(t>o{x) = 'Yh ^rnOX"^, ^o{y) = Y ^OnV'^, 

m>0 n>0 

m>l n>l 

It follows the functional equation Q{x,y)F{x,y) = G{x,y), which dehne the 
generating function F{x,y) = , where 

Q{x, y) = 1 - (7 + a)xy + ax'^y^, 

G{x, y) = -(7 a)xy{(l)Q{x) + '4)o{y) - yoo) - xyYn 

+y(t)i{x) + xifjiiy) + {(t)o{x) + ^poiy) - yoo). 

The polynomial Q{x, y) is the transform in ^ of the characteristic poly¬ 
nomial of the recurrence. 
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We write G{x,y) = K{x,y) — U{x,y), with 


K{x, y) = -(7 + a)xy{(j)o{x) + - > 00 ) - xyYn 


and 

U{x,y) = -y(t)i{x) -x^Jiiy) - ((j)o(x) +i’o(y) ->oo)- 

Then 

K(x,y) = F(x,y)G(x,y) + G(x,y), 

where K{x,y) is a polynomial in xy, like Q, whose coefficients are functions 
depending on the initial conditions. The function U{x,y) is affine, with 
coefficients depending on the initial conditions. 

Using the formal power series expansion 

00 

F{x, y) = G{x, y) ^ ((7 + a) - axyf x^y'", 

k=0 

we hnd the general term Ymn of the bivariate economic sequence. 
Particular solution If we use the following data 


(j)o{x) = loo 
^o(?/) = >00 
(j)l{x) = X 

^{y) = y 



Ymo = 0 for m > 1 

lo„ = 0 for n > 1 

111 = 1, Ymi = 0 for m > 2 

111 = 1, lin = 0 for n > 2, 


we obtain 

1 — (7 + a)xy + ax^y^ 

Denoting by ri,r 2 the real roots of the polynomial Q{x,y), it follows the 
bivariate sequence 




I f A B \ 
a ^ r”+V 


and Yrr, 


0 , for m ^ n, 


depending on the parameters a and 7 . 

Second computation variant Let us start with the univariate Samuelson- 
Hicks recurrence 

®n+2 — (o + 7)a„+i + attn = 0. 

Using a method similar to the previous one, we obtain the generating function 


F{x) 


do + («i - (a + j)ao)x 
1 — (a + j)x + ’ 
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where Oq, ai are fixed by the initial conditions. 

For the bivariate recurrence, the general solution, via characteristic equa¬ 
tion, is 

_ ^min{m,n} \min{m,n} 

mn — Cl Ai + C2A2 , 

where the coefficients ci, C 2 are fixed by initial conditions. Introducing the 
subsets Si = {m > 0, n > 0, m < n}, S 2 = {m > 0 , n > 0 , m > n} of N^, 
the generating function becomes 

F{x, y) = Y^ (ciXT + C2X^) (xyry^-^ + (ciA^ + c^A^) {xyTx^-^ 

Si S2 


E 

fc >0 


5 ^ (ciAr + {xyY 


m>0 


+ E 

k>l 


(cA? + c^a;) {xyY 


n>0 


x\ 


Using the univariate generating function, we find F{x,y) = , where 

Q{x, y) = 1- {a + 'j)xy + ax^y^, 

G{,x, y) = ^ [Yok + (hifc+i - [a + '^)Yok)xy] y^ 
k>0 

+ [Yko + {Yk+11 — (« + '))Yko)xy] x^. 


k>l 


5 Discrete multitime periodic coefficients 
Samuelson-Hicks model 


We have in mind the single-time case [S] and similar ingredients with those 
in Section 9. Instead of constant parameters 7 and a, we use variable pa¬ 
rameters 7 (f) and a(t), with 7 (f) G ( 0 , 1 ), a{t) > 0 , which are supposed to 
be T-diagonal-periodic. 

A discrete multitime periodic coefficients Samuelson-Hicks model can be 
either the first order diagonal recurrence system 


(Y(t + 1)A 

U(i+id 




subject to the (initial and compatibility) conditions flT^ . dH]), flT^ . or the 
second order homogeneous diagonal recurrence equation 


Y{t + 2) - ( 7 (f + 1) + a{t + l))Y{t + 1) + a{t + l)U(t) = 0, t G 
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subject to the (initial and compatibility) conditions flT7)) . flTS]) . flT^ . fEU]) . 

A prominent role in the analysis of the foregoing recurrences is played by 
the so-called Floquet multipliers of the second order homogeneous diagonal 
recurrence equation, if its coefficients are T-diagonal-periodic. We reconvert 
this recurrence into an equivalent diagonal recurrence system 

/ 0 1 A OVA 

\Z(t + l)) l,-a((+l) 7 (( + l) + a(( + l)t ' 

The matrix of this diagonal recurrence system 


A{t + 1) 


/ 0 1 \ 
y — (y.{t -h 1) -|- 1) -|- (y.{t -|- l)y 


is T-diagonal-periodic. A Floquet multiplier is nothing else than an eigen¬ 
value of the matrix 


T-l 

k=0 

This matrix is called monodromy matrix of the foregoing first order recurrence 
system. It is a 2 x 2- matrix with determinant 

A(f) = det D{t) = a{t + {T — ia{t)) ■ l) ■ ... ■ a(t + {1 — ia{t)) ■ l). 

Theorem 5. The Floquet multipliers of the second order homogeneous diago¬ 
nal recurrence equation, with diagonal-periodic coefficients, are the two roots 
of the quadratic equation — z{Tr D)-\-A = 0. These multipliers depend on 

For the constant coefficients recurrence equation, the foregoing equation 
is reduced to z^ — {^ + a.)z-\- a = 0 and the Floquet multipliers are constants. 

6 Conclusions 

This paper presents original results regarding the multivariate recurrence 
equations as continuation of [1]. Our approach to multivariate recurrence 
equations is advantageous for practical problems. The original results have a 
great potential to solve problems in various areas such as ecosystem dynam¬ 
ics, financial modeling, and economics. 

Though the techniques for deriving multivariate generating functions are 
sometimes paralleling the univariate theory, they achieve surprising depth in 
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many problems. Analytic methods for recovering coefficients of generating 
functions once the functions have been derived have, however, been sorely 
lacking. 

The authors lay no claims to the paper’s being a complete presentation 
of all current methods for investigation the linear multivariate recurrences 
with periodic coefficients. The material presented here is a reflection of our 
scientihc interests regarding Floquet theory. 

We hope to recruit further researchers into multivariate recurrences, which 
still have many interesting challenges to offer, and this explains the rather 
comprehensive nature of the paper. 
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